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Abstract. Yor's generalized meander is a temporally inhomogeneous modification of the 2(j/+l)-dimensional 
Bessel process with v > —1, in which the inhomogeneity is indexed by n € [0, 2(y + 1)). We introduce the 
non-colliding particle systems of the generalized meanders and prove that they are Pfaffian processes, in the 
sense that any multitime correlation function is given by a Pfaffian. In the infinite particle limit, we show 
that the elements of matrix kernels of the obtained infinite Pfaffian processes are generally expressed by 
the Riemann-Liouville differintegrals of functions comprising the Bessel functions J v used in the fractional 
calculus, where orders of diffcrintegration are determined by v — k. As special cases of the two parameters 
(u, k), the present infinite systems include the quaternion determinantal processes studied by Forrester, 
Nagao and Honner and by Nagao, which exhibit the temporal transitions between the universality classes of 
random matrix theory. 
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1 Introduction 

The random matrix (RM) theory was introduced originally as an approximation theory of statistics of nuclear 
energy levels j^O]. It should be noted that at the same time as the standard theory was established for three 
ensembles called the Gaussian unitary orthogonal, and symplectic ensembles (GUE, GOE, GSE) JT], Dyson 
proposed to study such stochastic processes of interacting particles that the eigenvalue statistics of RMs are 
realized in distribution of particle positions on R |10| . Dyson's Brownian motion model is a one-parameter 
family of A-particle systems, Z (/3) (i) = (z[ p \t), zf \t), • • • , Z$\t)), described by the stochastic differential 
equations 

7^ p) (i\ 7 yp) (i\ 
l<j<N,j=^i W ~~ A j W 



dZ^\t) = dB i {t) + ^ ^)7T 75*77*' te[0,oo),l<i<N, (1.1) 

2 . . i rr? . , . 7 (P> - 7S P) ^ 



where Bi{t),i = 1, 2, • • • , N are independent standard Brownian motions and the parameter (3 equals 2, 1 
and 4 for GUE, GOE and GSE, respectively. Due to the strong repulsive forces, which are long-ranged 
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and act between any pair of particles, intersections of particle trajectories are prohibited for (3 > 1 43 
(see also [J]). In this one-parameter family, the (3 = 2 case (i.e. the GUE case) is the simplest and the 
most-understood, since its equivalence with the N particle systems of Brownian motions conditioned never 
to collide with each other can be proved |16|. 

The standard (Wigner-Dyson) theory has been extended by adding three chiral versions of RM ensembles 
in the particle physics of QCD |52l 1511 ITiH I4fi| . and by introducing the four additional ensembles so-called 
the Bogoliubov-de Gennes classes in the mesoscopic physics Q]|2]- Here we note that the chiral ensembles 
have a parameter v G {0, 1, 2, • • • } in addition to (3. In these totally ten ensembles ^ |^] [2], chiral GUE 
(chGUE), class C and class D can be regarded as natural extensions of the GUE, in the sense that these 
eigenvalue statistics are also realized in appropriate non-colliding systems of stochastic particle systems: 
Konig and O'Connell showed that the chGUE with the parameter v S {0, 1, 2, ■ • • } corresponds to the non- 
colliding systems of 2(v + l)-dimensional squared Bessel processes |28| . The present authors clarified that 
the eigenvalue statistics in the classes C and D are realized by the non-colliding systems of the Brownian 
motions with an absorbing wall at the origin and of the Brownian motions reflecting at the origin |271 126| . 
Since the absorbing and reflecting Brownian motions are directly related with the three-dimensional and 
one-dimensional Bessel processes, respectively (see, for example, @U), the stochastic differential equations 
of these non-colliding particle systems are generally given by 

2^ + 1 1 x ^ f 1 1 \ 

2 z?\t) + x < 1 zt\t) - zf\t) + z\ v \t) + zf\t) J J ' 

tG[0,oo),l<i<JV, (1.2) 

with reflecting barrier condition at the origin in case v = —1/2. Therefore, the difference of (non-standard) 
RM ensembles can be attributed to the difference of dimensionality of the Bessel processes, whose non- 
colliding sets realize the statistics of the RM ensembles Here we remind that the d-dimensional Bessel 
process is defined as the process of the radial coordinate (the modulus) of a Brownian motion in M. d . To realize 
other 10 — 4 = 6 RM ensembles by conditioned stochastic processes may be much more difficult (see |49| ) . 
but we demonstrated that, if we consider appropriate non-colliding systems of temporally inhomogeneous 
processes defined only in a finite time-interval [0,T], we can observe the transitions of distributions into the 
6 distributions as the time t approaches the final time T [2S1 I2EI • The interesting fact is that the processes 
that can be used instead of the Bessel processes 111 .21) should have one more parameter k in addition to v. 
This two-parameter family of temporally inhomogeneous processes indexed by (y, k), v > — 1, k . G [0, 2(y+X)) 
is equivalent with the family of processes already studied by Yor. He called them the generalized meanders 

m 

From the view-point of random matrix theory, studying time-development of stochastic systems by calcu- 
lating, for example, the multitime correlation functions corresponds to considering multi-matrix models. In 
particular, the temporally inhomogeneous processes will be identified with such matrix models that matrices 
with different symmetries are coupled in a chain |231 1241 1341 122] . Determination of all multitime correlation 
functions of systems, which allows us to determine scaling limits associated with the infinity limit of matrix 
sizes (i.e. the infinite-particle limit) is one of the main topics of the modern theory of RM [1231 ■ The finite and 
infinite particle systems showing the orthogonal-unitary and symplectic- unitary transitions, and transitions 
between class C to class CI were studied and multitime correlation functions were determined by Forrester, 
Nagao and Honner (FNH) To] , and by Nagao [32], respectively. The system in the Laguerre ensemble with 
(3 = 1 initial condition reported in the former paper can be regarded as the v = k G {0, 1, 2, ■ • • } case of 
the non-colliding system of the generalized meanders and the system reported in the latter paper as the 
(v, k) = (1/2, 1) case. 

If we think about the system of generalized meanders apart from the RM theory, however, we can consider 
the parameters v and k as real numbers, and not necessarily integers nor half-integers. In the present paper, 
we calculate the multitime correlation functions of non-colliding systems of (squared) generalized meanders 
for arbitrary values of parameters, provided they satisfy the condition v > —1, k G [0, 2(v + 1)) so that the 
systems are not collapsed. We first define the N particle systems in a finite time-interval [0, T] and take the 
N = T — > oo limit to construct the two-parameter family of infinite particle systems. We prove that the 



dZ\ v \t) = dBi(t) + 
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multitime characteristic functions is given by a Fredholm Pfaffian JJ| and thus any multitime correlation 
function is given by a Pfafhan. Similarly to the results by FNH and Nagao [221 and their temporally- 
homogeneous version (the determinantal process with the extended Bessel kernel IH3); the elements of the 
matrix kernels of Pfafhans are expressed using the Bessel functions, but we clarify the fact that they are 
generally given by the Riemann-Liouville differintegrals of the functions comprising the Bessel functions, 
which are used in fractional calculus (see, for example, |3fil 1441 l3"%] L This structure will explain the origin 
of the multiple integral expressions for the elements of the matrix kernels reported by FNH J^l and Nagao 

The paper is organized as follows. In Section 2, the definitions of the generalized meanders of Yor and 
their non-colliding systems are given and the Riemann-Liouville differintegrals of the Bessel functions with 
appropriate factors are introduced. The main theorem for the infinite particle limit (Theorem 12. 1J) is then 
given. It is demonstrated that, if we take a further limit in the system of Theorem 12.11 we will obtain 
the temporally homogeneous system of infinite number of particles, which is a determinantal process with 
the extended Bessel kernel studied in [SU| (see also [37]). Using the properties of the Riemann-Liouville 
differintegrals, we show that Theorem 12.11 includes the results by FNH [IS] and Nagao [321 as special cases. 
Section 3 is devoted to prove that for any finite number of particles N, the present system is a Pfaffian process 
fThcorcm l3.1|) , in the sense that any multitime correlation function is given by a Pfaffian 0Oj • These Pfaffian 
processes may be regarded as the continuous space-time version of the Pfaffian point processes and Pfaffian 
Schur processes studied by Borodin and Rains Soshnikov used the term Pfaffian ensembles in |Sl l47ll4*5] , 
See also [351 HOI ESI H31 ITT) in the context of study of nonequilibrium phenomena in the polynuclear growth 
models, and |35M14| in that of shape fluctuations of crystal facets. The processes studied in ^3E2] are also 
Pfaffian processes, since the 'quaternion determinantal expressions' of correlation functions, introduced and 
developed by Dyson, Mehta, Forrester, and Nagao |12l 1291 l3T)l 1331 131j . are readily transformed to Pfaffian 
expressions. The method of skew-orthogonal functions associated with the Laguerre polynomials are used 
in Section 4 in order to perform matrix inversion and give explicit expressions for the elements of matrix 
kernels of Pfaffians. Asymptotics in T = N — > oo are studied in Section 5. Appendices are given to show 
proofs of formulae and lemmas used in the text. 

At the end of this introduction, we would like to refer to the papers [IHIIE]> which reported the further 
extensions of RM theory in physics and the representation theory. We hope that the present paper will 
demonstrate the fruitfulness of developing the probability theory of interacting infinite particle systems in 
connection with the extensive study of (mufti-) matrix models in the RM theory. 



2 Definition of Processes and Results 

2.1 Non-colliding systems of generalized meanders 

Let Z and R be the sets of integers and real numbers, respectively, and set N = {1,2,...}, No = NU {0}, 
Z_ = Z \ N , and R+ = {x € R : x > 0}. Let T(c),c € R \ (Z_ U {0}), be the Gamma function: 
r(c) = / °° dy e-yy'- 1 for c> 0, and T(c) = T(c + [-c] + l)/{c(c + 1) • • • (c + [-c])} for c € (-oo,0) \ Z_, 
where [c] is the largest integer that is less than or equal to the real number c. For t > 0, x,y £ M. + and 
v > — 1 we denote by G^\t] y\x) the transition probability density of a %y + l)-dimensional Bessel process 

sum, 

GW(t;y\x) = ^e-^+y^lJ^), X >0,yeR + , 
x v t \ t J 

G { "\t-y\Q) = f V+ \ ^ e- y2/2t , y £ R+, 

where I v (z) is the modified Bessel function : I v (z) = E^o( z / 2 ) 2n+ "/{ r ( fl + 1 ) T ( V + n + For T > °» 
k £ [0, 2(u+ 1)), we put 

hp K) (t,x)= dyG^(T~t;y\x)y~ K 7 x € R+, t e [0, T], 
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and 

Gp K \s,x;t,y) = ^ 77 ±- -G«(t- s;y\x)hp K \t,y), x>0 7 yeR +} (2.1) 

h T ' (s,x) 

Gp K \0,0;t,y) = vf T [\ + 1} ... (27y/ 2 G<"> (t; y|0)4^ ft, y), j, e R + , (2.2) 
1 [y + 1 — k/2) 

for < s < t < T. This transition probability density G^''^ (s, x; t, y) defines the temporally inhomogeneous 
process in a finite time-interval [0, T], which is called a generalized meander. In particular, when v = 1/2 
and K = 1, it is identified with the process called a Brownian meander (see Chapter 3 in Yor |53|). 

Now we consider the iV-particle system of generalized meanders conditioned that they never collide in a 
finite time-interval [0, T]. Let 



= jx = (zi, #2, • • • ) xjv) G :0<xi <X2 < ••• < Xn\- 



According to the determinantal formula of Karlin and McGregor the transition probability density is 
given as 

g K N ' T ; (s,x;t,y) = -^y- , < s < t < T, x,y£l +< , (2.3) 

■^n't ( T ~ s > x ) 

where 

/fe^s.xj^y) = det Gp K ' (s,Xj,t,y k ) , A/T'^ (i,x) = / dyf^'(T - t,x;T, y). 
l<J>fc<JV L J Jr^ 

Since h^'°\t,x) = 1, G^ (s,x;t,y) = G^ v \t — s;y\x) and thus /j^'^ is temporally homogeneous and 
independent of T, we will write (t — s;y|x) for f^'j) (s,x;f, y). Moreover, note that 

i^( s ,x;i,y) = -TTTT^; — r /^ft- S ;y|x)^ K) (t,y), 

/ly' (s.x) 



where hj, (i, x) = Yif=i hj! (t,Xj) and /^'^(T, x) = YljLi x j " ■ Then l|2.3|) can be written as 

1 

a71T' k) (t- s ,x) 



x;t,y) = -^y^ :/#>(* - s; y\x)X#' K) (T - t, y), (2.4) 



where 

AT 



A/lrV,x) = / rfy /lr } (t;y|x) (2.5) 

In our previous paper |2fi| it was shown that, taking the limit x — > = (0, 0, . . . , 0) at the initial time s = 0, 
l|2.4|) becomes 



JV 



(0, 0; t, y) = C^ T (i) J] (t, |0) J[ (y 2 k - 2/f) A#' K) (T - t, y) 

i=l l<j<k<N 



(2.6) 



for v > —1 and k G [0, 2(f + 1)), where 



_ i r i> + i)r(i/2) 

U N,TK t ) - oJVW-/t-l)/2 11 



2 w(jv- k -i)/2 JUL r (j/2) r ((j + 1 + 2v - k)/2) ' 
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The TV-particle system of non-colliding generalized meanders all starting from the origin at time is defined 
by the transition probability density g^'ft given above and it will be denoted by X(t) € M/^ < ,t 6 [0, T] in 
the present paper. It makes a two-parameter family of temporally inhomogeneous processes parameterized 
by v > -1 and re e [0, 2(i/ + 1)). 

We denote by X the space of countable subsets £ of K satisfying n if) < oo for any compact subset X. 
For x = (xi,^, ■ ■ ■ ,Xn) S Ufci R£ : we denote {x t }f =1 € X simply by {x}. Then 3^(i) = {X(*)},t € [0,T], 
is the diffusion process on the set 3£ with transition density function (s, £; t, 77), < s < t < T: 

f ff^(a,x;t,y), if s > 0, jj£ = ftf? = iV, 

flfer(*.f;*.»7) = S ^(O.Oji.y), if s = 0, £ = {0}, (jry = TV, 
[ 0, otherwise, 

where x and y are the elements of R+< with £ = {x}, r\ = {y}. 

For the given time interval [0, T], we consider the M intermediate times < t\ < ti < ■ ■ ■ < tM < T. 
For convenience, we set t = 0, t M +i = T. For x< m ) e R N , 1 < m < M + 1, and TV' = 1,2,...,JV, we 

put x^ = f a;^ 71 ^, X2™\ • • • , and £^ = {x^}. Then the multitime transition density function of the 

process S^(t) is given by 

M 

where we assume £^ = {0}. For a sequence {-^V m }m=i of positive integers less than or equal to N, we define 
the (iVi, -/V2, . . . , A r M+i)-niultitime correlation function by 



X /. t JVi , e N 2 . .+ C N M +: 
PN,T l r l'Sl : r 2,? 2 1 • • • ! IM+l, SjVf+1 

M+1 . N 



I II iN \ V II ^ fl ^(Q,{0} ;il) ef;...;iM + i,eM + i)- (2- 

• — =1 ^ m '' J=N m + l 



Associated with the generalized meander <|2.1II , (|2.2|) , we consider a temporally inhomogeneous diffusion 
process with transition probability density 



pp K \s,x;t,y) = G^' k) (s, yfx : t, y/y) x -y 1 / 2 , x>0,yeM.+ , 

t G [0,T], and call it a squared generalized meander. The ./V-particle system of non-colliding squared gener- 
alized meanders Y(i),t £ [0, T], is then defined by 

Y(t)= (^(t^.Xaft) 2 ,...,^*) 3 ), te [0,T]. 

The correlation function p)v T of = {Y(i)} is obtained from H2.8(l through the relation 

Y A aJVi . . /-N 2 . . . /-N M + i 

PN,T I r l > SI ! r 2,t, 2 , • • • )tM+lj S>M+1 



JV,T (il^f 1 ;^,^ 2 ! ■ ■ ■ ^M+l^M+l 1 ) II II (mj> ( 2 - 9 ) 

m=l j=l 



where = {xg}, = {yg> with = 0/j m) , 1 < j < N m , 1 < m < M + 1. 
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2.2 Riemann-Liouville differintegrals of Bessel functions 

We consider the following left and right Riemann-Liouville differintegrals for integrable functions / on R + , 

1 / d x " rx 



nif(x) = jl^j J o (x~yr- c - 1 f(y)dy, (2.10) 

where cgl and n = [c + 1]+ with the notation x+ = max{i, 0}. It is easy to confirm that, if c £ No, both 
of them are reduced to the ordinary multiple derivative, 

and, if c e Z_, they are equal to the multiple integrals, 

o^>xf( x ) = / dj/M_i / dj/|c|-2 ' ' ' / ^2/1 / dy Q f(y ), 
Jo Jo Jo Jo 

/•CfO poo />oo poo 

x® c oof{x) = / dy|c|-i / dy\ c \_ 2 ■■■ dyi dy f(y ). 
Jx Jy\c\-i Jyi Jyi 

For c € (-co, 0) \ Z_ (I2.1U|) and l|2.11[l define fractional integrals, and for c € K + \ No fractional differentials. 
The Riemann-Liouville differintegrals are most often used in the fractional calculus (see, for example, 36 

mmi). 

Let J v (z) be the Bessel functions: J u (z) = Y.7=o(~ 1 ) 1 ( z 1 2 ) 2l+U I i T ( v + 1 + l)t\}. We define functions J u 
and J v as 

= (^r /2 J,(2^)e 2 ^ = e 2 ^£ tlfe^, (2 . 12) 

OO 



1(9, n, x, s) = (9 VX )-^M2^)e^ = e^> £ -XM_ (2 . 13) 



We will use the following abbreviations for the Riemann-Liouville differintegrals of order c £ K of J y and J„, 

jW(6,ri,x,s) = T> c n J v (6,r),x,s), 6,n>0,seR, (2.14) 

j( c >(M,x, S ) = „D^J„(#, tj.i.s), 9, V >0,s<0. (2.15) 

We note that, if c e IR \ No, Ji can be expanded as 

Jj > °\9,ri,x,8) = w - -J2 u V \ Ji n) (9,V,x,s), 9, v >0,seR. (2.16) 
1 — c * — ' n (n — c 

It is also noted that, since J v (9, rj, x, s) — > exponentially fast as rj — ► oo, if s# < 0, 

jM(0,77,x, S )-- -/ de^-Tj)"- - 1 ^^^,^*), M>0,s<0, (2.17) 

r(n - c) J v 

where n — [c + 1] + . 
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2.3 Results 

We put 



a = a(f, k) = v — — , fa = b{v, k) = v — k, 



(2.18) 



and introduce functions T>(s, x; t, y), l(s, x; t, y), S(s, x; t, y) and S(s, x; t,y), x, y G R + , s, t < 0, 

1 r 1 r~ 

v(s,x;t, y ) = d ee 1 - K [4- b - 1 \e,i,x,-s)4- b \e,i,y,-t) 

-ji- h \6,l,x,-s)Jl- b - l \6,l,y,-t) , 

/OO /*00 
dfl r- 1 [ J d£ eZ b+1) (0, e, a:, a) # b+1) (0, !,»,*) 

/■OO 

-# 6+1 >(M,*,*) / d££°# b+1 >(0,£,j/,i) 



S(s,x;t,y) 



1 /a; 



re/2 „i 



d0 



^ 6+1 )(e,l,a; ) *)4- 6 - 1 )(6l,l,y,-t) 



and 



2 Vy/ Jo 

/OO 

S(s, x; t, y) = S(s, x; t, y) - l (s<t) ^-J Q{s, x; t, y), 
where is the indicator function: = 1 if uj is satished and =0 otherwise, and 

/•OO 

a:; t,i/)= / dQ J v (2^)J v (2^)e 2 ^ e . 
Jo 

For an integer N and a skew-symmetric 2N x 2iV matrix A = (ay), the Pfaffian is defined as 

Pf(A) = Pfl<j<j<2Ar(ajj) = ^2 s S n ( (T ) a CT(l)(T(2) a CT(3) CT (4) ■ • • a C r(2JV-l)o-(2Ar) ) 



(2.19) 



(2.20) 



(2.21) 



(2.22) 



where the summation is extended over all permutations a of (1,2, . . . , 2N) with restriction <r(2k — 1) < 
a(2k),k = 1,2,..., N. We put 

Sjv(«) -{^l(T A r + s),y 2 (T A r + S ),...,y A r(T A r + s)}, St [-Tjv,0), 

and S^(s) = {0}, s G (— oo, —TV). Then we can state the main theorem in the present paper. 

Theorem 2.1 Le£ Tjv = JV. T/ien i/ie process E^(s),s G (— oo, 0) converges to the process E^(s),s G 
(—oo,0), as N — > oo, m £/ie sense of finite dimensional distributions, whose correlation functions p Y are 
given by 

p* {y g}; S2 , {y<g}; ; {yg?}) = Pf [-4 (y<g , y«, . . . , y^ , 

/or any M > 1, any sequence {N m }^ =1 of positive integers, and any strictly increasing sequence {s m }m=i 
of nonpositive numbers with s M +i = 0, where A (y$, y^, . . . , y^) is the 2^" =1 iV m x 2^" =1 iV m 
skew-symmetric matrix defined by 



(A m - n (v (m) v {n) )) 



<j<AT„,l<m,n<M 



with 2x2 matrices A m,n (x,y) ; 



A m > n (x,y) 



In the infinite-particle system defined by Theorem l2.ll we can take the further limit: 

s m — > — oo with the time differences s n — s m fixed, 1 < m, n < M. 

In this limit, D(s m ,x;s n ,y)T(s m ,x;s n ,y) — * 0, 1 < to, n < M, as we show in Appendix El Therefore, we 
can replace I? and X by zeros in the matrices. Then the Pfaffian is reduced to an ordinary determinant of 

the Z%±}N m x Y^N m matrix, A (yg.yg, . . . ,y^) = (a-"^,^))^ 

with the elements 



l<i<AT m ,l<j<Af„,l<m,n<M 



(m) _ (ri) 



where 



2( s -t)e 



( f 1 

/ dO J^(2v / fe)J„(2 v /%)e 2 
S( 8 ,x;t,y)= | J vPVZ)y/yJU2Vv) ~ M2^y)^J' v {2^) 



x-y 



rOC ____ 

- / d0 J v (2\f0x)J v (2^/6y)e 
Ji 



2(s-t)0 



if s > t, 



if s = <, 



if s < t, 



with J^, = dJ u {z) / dz. Hence, in this limit we obtain a temporally homogeneous system of infinite number 
of particles, whose correlation functions are given by 



P Y {y$}; s 2 , {y«}; . . . ; {y<£}) = A (y« , y<g, . . . , y^ } ) - 



(2.23) 



Remark 1. Forrester, Nagao, and Honner ^5] studied the orthogonal- unitary and symplectic-unitary 
universality transitions in random matrix theory by giving the quaternion determinantal expressions of (two- 
time) correlation functions for parametric RM models. One of their results for the 'Laguerre ensemble with 
0=1 initial condition', which shows the orthogonal-unitary transition, can be reproduced from Theorem 
I2.1l bv setting 

(i) k = v -i=^> a= 2' k = ^' where i/ £ No. 
This fact may be readily seen, if we notice that by definition 



4- x \9,l,x,a) 



dr, (6r]x)" /2 J42y/(hjx~)e 2seT1 



^x^ 2 / duu^ 2 J„(2V^x~)e 2su . 



Remark 2. Nagao's result on the multitime correlation functions for vicious random walk with a wall 32 
can be regarded as the special case of Theorem 12 .11 in which 

(ii) v= \ K = 1 a = ' b = ~\- 



S 



This fact can be confirmed by noting that, by definition Ij2.14|l with J\/2{6, 0; x, s) = 0, 

(te) 1 /* r 1 



^1/2 1; 5 ) 



d7 7 (l-7 ? )- 1 / 2 r, 1 / 4 J 1/2 (2V^)< 



2s0r, 



(to) 1 / 4 r 1 



d»> (1 - vr 1/2 -^ {?1 1/4 J 1/2 (2^)e 2 ^} , 



by (HZ), 

and, by definition Q2.15J1 . 

■?(l/2) 



* (C ~ {r 1/4 J 1/2 (2^)e 2se «} , a < 0, 



-1/4 



rf7/(77 — l)- 1/2 T ] - 1 / i J 1/2 {2 v ^)e 2se \ .s ( I. 



In this case, the system shows the transition between the class C and class CI of the Bogoliubov-de Gennes 
universality classes of nonstandard RM theory |321 1271 . 

Remark 3. From the results for finite non-colliding processes [2H|j we expect that, when 

v — 1 



(iii) k = v + 1 



a = 



b = —1, where v e No, 



the present infinite particle system will show the transition from the chiral GUE to the chiral GOE of the 
universality classes and when 



(iv) v = ~\ K = 



o = b = -- 



1 



that from the class D to the 'real-component version' of class D of the Bogoliubov-de Gennes universality 
classes |2~fi] . 

Remark 4. Following the argument given in |39l I2()j . tightness in time can be proved and transition 
phenomena observed in the limit sm —> may be generally discussed, which will be reported elsewhere. 

Remark 5. The homogeneous system H2.2rt|l was studied in |50l 137) . 

3 Correlation Functions Given by Pfaffians 

3.1 The multitime transition density 

If we put 

G^ K \t,y\x)=G^(t,y\x)x (|)~ K , x>0,yel + , 
&^(t,y\0) = G^(t,y\0) x y~«, y G R+, 
the multitime transition density (|2.7|) with to = 0, tjvf+i = T, and £o — {0} is written as 

fl^ (O, {0}; h, {xW}; • • • ; t M+ u {x^ +1) } 



n {(4 i3 ) 2 -H 13 ) 2 } n ^ M+1 



l<j<k<N 



l<j<k<N 



N 



M 



TTg^^i.^-Io) TT dct 

11 3 LL l<j,k<N 



3=1 



G^)(t m+ l-*m,4 m+1) |4 m) ) 



where Ij2.4|l and Ij2.6|l with (|2.5[) are used. 

Through the relation H2.9|l . the multitime transition density for the process {Y(t)},t £ [0, T], denoted 
by pt'm is then written as 



where 



;f(0,{0};t 1 ,{yW};...;W,{y (M+1) }) 

<^ T (OMy (1 W(My (M+1) )) 

JV M 

n^^.^io) II ; .<i<;< v 



k=l 



^'(Wi-U^ljf" 1 ) 



^Ar(y) = n {yj-Vi), ye 

l<i<j<JV 

f^ K >(i,y|0) = G^(t,^|0)xiy- 1 /2 



siV 



(3.1) 



-(x+2/)/{2(i- S )} 



(f) 



6/2 



2(t-s) Vx 

Expectations related to the process {Y(ti)}, {Y(t2)}, • ■ ■ , {Y(ijv/+i)} are denote by T 

M+l r . M+l 



x > 0, ye K+. 

Y 



N.T 



f({Y(t 1 )},{Y(t 2 )},...,{Y(t M+ i)}) 



1 



n^y 



(m) 



rn — 1 



(3.2) 



x/({y (i n,{y (2J },---,{y (M+1) }^ (3-3) 

3.2 Predholm PfafRan representation of characteristic function and PfafRan pro- 
cess 

For simplicity of expressions, we assume from now on that the number of particles N is even. The references 
jS3E2 will be useful to give necessary modifications to the following expressions in the case that N is odd. 
Let Cq(M) be the set of all continuous real functions with compact supports. For f = (/i, /z, • • • , Jm+i) £ 
Co(M) +1 , and 6 = (81,82,- ■■ ,&m+i) £ M M+1 , the multitime characteristic function is defined for the 
process {Y(t)},t £ [0,T] as 



{M+l JV ~j " 

m— 1 i m -l ) _ 



(3.4) 



Let x m (x) = e ^r^-S m f m (x) — 1, 1 < m < JWf + 1- Then by the definition of multitime correlation function 
lf2~3|) with ffity . we have 



TV JV 



JV M+l 



JVi=0 JV 2 =0 JVm+1=0 m=l 
M+l JV, 

X 

m=l iim) = i 



(M+l) 
JV A/ + i 



II II Xm^yl^PNp (*i.{y^ ) }; i 2,{y^};---;iM+i,{yi^+i ) } 



(3.5) 
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that is, the multitime characteristic function is a generating function of multitime correlation functions T . 
We consider a vector space V with the orthonormal basis \ \m,x) > , which satisfies 

I J l<m<Af+l,x€K+ 

(m,x\n,y) = 5 mn S(x-y), m,n = 1, 2, • • • , M + 1, x, y e R+, (3.6) 

where <5 m „ and <5(x — y) denote Kronecker's delta and Dirac's 5-measure, respectively. We introduce the 
operators J,p,p + ,p_ and \ acting on V as follows 

(m,x\J\n,y) = l (m= „ =M+1) sgn(y - x), (3.7) 

(m,x\p\n,y) = l( m <n)P^ v ' K \t n -t m ,y\x) + l {m>n) p^ v ' K \t m -t n ,x\y) 

+ l {m=n) 5(x-y), (3.8) 

(m,x\p+\n,y) = l( m <„)P ( ' / ' K) (t„ - t m , y\x) = (n, y\p-\m, x), (3.9) 

(m,x\x\n,y) = Xm(x)S mn S(x ~ y), (3.10) 



and we will use the convention 

M+l 



(m,x\A\n,y){n,y\B\£,z) = ^2 / dy A(m,x;n,y)B(n, 

n=l 



j/;f,z) = (ro,a;|AB|*,z> 



for operators A and i? with (m, x|A|n, y) = A(m, x; n, y) and (m, x\B\n, y) = B(m, x; n, y). 

Let Mi(x) be an arbitrary polynomial of x with degree i in the form Mi{x) = biX 1 + ■ ■ ■ with a constant 
bi 7^ for i e Nq. Since the product of differences /i/v(x) is equal to the Vandcrmonde determinant, we have 



w fe=l 



m^htt^-i? det ^-i(^) • ( 3 - n ) 

I A A J l<i,j<N L 



Then we consider the set of linearly independent vectors ; i e n| in V defined by 

|i) = \m, x)(m, x\i), 

where 



(m,x\i) = (i\m,x)= / dyMi-id/^^^.ylO^"^^-*!^!!/), (3.12) 
i e N, m = 1, 2, . . . , M + 1, a; € R+. We will use the convention 

oo 

<i|A|j)0-|S|m >a: > =^A ij sj ro) ( a; ) - (*|AoB|m,x), 



are not assumed 



for Aij = (i\A\j) and B^ n \x) — (j\B\m, x). It should be noted that the vecors ||«) ; i e n| 

to be mutually orthogonal. By these vectors, however, any operator A on V may have a semi-infinite matrix 

representation A = ( (i\A\j) ) . If the matrix A representing an operator A is invertible, we define the 

V / i,jeN 

operator A A so that its matrix representation is the inverse of A; 

((i\A*\j)) =A-\ (3.13) 

that is, (i\A\j)(j\A A \k) = (i\AoA A \k) = Sit, i.fce N. 
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Let Vn be a linear operator projecting 5pan||i) ; i G Nj to its iV-dimensional subspace Span^ 
1,2, . . . ,7V | such that 

{i\m,x}, ]fl<i<N, 



i) ; i 



(i\V N \m,x) = (m,x\V N \i} 



otherwise. 



We will use the abbreviation An = VnAVn for an operator A. If the N x N matrix defined by An = 
({i\A N \j))i<i,j< N is invertible, then (A N ) A is defined so that ( (i\(A N ) A \j) J = (An)^ 1 , and (i\(A N ) A \j) = 

\ / I <i,j< .v 

0, if i > iV+ 1 or j > N + 1. 

As shown in^ we can prove that 



*y T (f . ) I = Det I 2 <5 m „<5(x - y) + 



S* m '"(x,y) I m - n (x,y) 
D m ' n (x,y) S n ' m (y,x) 



Xn(y) 



where Det denotes the Fredholm determinant. Here I 2 is the unit matrix with size 2, 

D m ' n (x,y) = -(m,x\o(j N )*o\n,y), 
S m > n (x,y) = (m,x\pJo(j N ) A o\n,y), 
I m > n {x,y) = -(m,x\pJo(j N ) A oJp\n,y), 

S m < n (x,y) = S m ' n (x,y)~(m,x\p + \n,y) 
I m ' n (x,y) = I m > n (x,y) + (m,x\pjp\n,y). 

It implies that the multitime characteristic function is given by the Fredholm Pfaffian |4Uj . 



and 



*£ iT (f; 0) - PF(j 2< W(* - y) + ^Jx)A m ^{x, y)V)Uy)) , 



where J 2 



1 

-1 



and 



J 2 



S" m '"(a;,y) I m ' n (x,y) 

D m ' n (x,y) S n ' m {y,x) 

D m ' n (x,y) S n ' m (y,x) 

-S m ' n (x,y) -I m > n (x,y) 



It is defined by 



PF [J 2 S mn 5(x -y) + VxMA m - n (x iy )y/^j 

N N N M+l 



*!2 



EE- E 

Ar 1= 0Af 2 =0 iVj, +1 =0 m=l 
M+l N m 

n n x m (»£0 K (^A---^) 



dyiv M+1) 

17 ^ M + l 



m— 1 i(m) = i 



(3.14) 



(3.15) 

(3.16) 
(3.17) 



(3.18) 



(3.19) 



where A (y%\ , y^ , • • • ,yjy£^) denotes the 2 X)m=i N ™ x 2 X^^=i -^m skew-symmetric matrices con- 
structed from H3.18JI as 



L(ySg.ySr--,yl£t?) = (^ m) ,i ( j" ) )) 



l<i<AT m ,l<j<Af„,l<m,n<M+l 
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for N m = 1, 2, • • • , N, 1 < to < M + 1. Comparison of 1|3.5|1 and (j3.17|l with H3.19J1 immediately gives the 
following statement. 

Theorem 3.1 The N -particle non- colliding system of squared generalized meanders Y(t),t S [0, T] is a 
Pfaffian process, in the sense that any multitime correlation function is given by a Pfaffian 

pIt (it, {y£>}; h, {yg}; . . . ; W, {y^ 3 }) = Pf ^(yg.yg, • • • , C)) • 



4 Skew- Orthogonal Functions and Matrix Inversion 
4.1 Skew-symmetric inner products 

Consider the N x N skew-symmetric matrix Aq = ((^4o)«j)i<i.j<A r with 

(A )ij = (i\JN\j) = (i\rn,x)(m,x\J\n,y)(n,y\j), i,j = 1,2, ... ,N. 



(4.1) 



In order to clarify the fact that each element (Aq)^ is a functional of the polynomials Mj_i(:r) and Mj_±(x) 
through 1)3. 12[) . we introduce the skew-symmetric inner product 

/•oo /*oo 

(f,g}= dx dy F(x,y)p<^(t 1 ,x\0)^Ht uy \0)f(x)g(y), 



where 

Fix v) - Hdw l W dz $ y ' K \T-*uz\x) ^\T-t u w\x) 
Then we have the expression 

(A )y = (Mi_i,Mj_i), '../ 1.2.---..V. 
We now rewrite the skew-symmetric inner product i|4.2|) by using the simpler one 

(/,<?>* = -<<?,/>* 

rfwe^^w / dze- z/2 z a J/(z) ff ( W )-/H 5 (z) 



(4.2) 
(4.3) 
(4.4) 



(4.5) 



which we call the elementary skew-symmetric inner product. Remind that p( u ' K *> is given by 13.21) us- 
ing the modified Bessel function. We will expand it in terms of the Laguerre polynomials, L"(x) = 
(x~ a e x / jl)(d/ dxY (e~ x xi +a ) , a S M, j G No, using the formula 



E 

i=o 



r(j + l)Ly(a;)Ly(yV 
r(j + 1 + 1/) 



< 1,1/ > -1. 



(4.6) 



(See the corresponding calculation for the non-colliding Brownian particles in [221, where the heat kernel 
was expanded in terms of the Hcrmite polynomials.) For this purpose, it is useful to introduce the variables 



tji{2T t n ) 

f : 



2T — t 

Xn = " , n = l,2,...,M+l, 



since we can see that 



2(t n -t m ) y l-Xn/Xm / V c m? 

i+M 



b/2 



x exp 



1 - Xn/Xr, 



2T J 



1 - Xn/Xr 



2T ' ' ; 
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and, if we apply the formula (|4.6(1 with r = x n /Xm, ^ = £ and y = 77, it is written as 



^ in 

00 

E 

j=0 



c-°- 1 e /2 (c„77) a exp 



r(j + 1) 
ro' + i + i/) Vx 



— — £ — d — — ) 

2T [ 2T )T] 



(4.7) 



That is, c n and give the spatial scale of spread of TV particles and the proper temporal factor at time 
t n , respectively. (See equation (17) and explanation below it in [12], where the variable c„ was determined 
by showing that the one-particle density obeys Wigner's semicircle law scaled by c n for the non-colliding 
Brownian particles.) In particular, for n = M + 1 we have 



^\T-t m ,T V \c m O 



t: 



J 1 k/2+1 



-o— ltft/2 a 



r ; V exp 



xe-«/ 2 e- 



" /2 E 

J=0 



ECz + 1) 

T(j + l + v 



T 2 



(4.8) 



since cm+i — T and xm+i = 1- Then we obtain the relation 



Cl 



Cl 



r> + i) 2 y 



oo oo 



cfa / dye- x e-yx»y v f{x)g{y) 



o 



j=0 fc=0 



r(i + i) _ x „ r(fc + i) 



rtf + i + i/) J T(fc + i + z/) 



(4.9) 



4.2 Skew-orthogonal polynomials 

For a£l and n £ Z we define 

r(n + a + 1) 



n + a 
n 



r(n + l)r(a + 1) 
(_l)»r(-a) 

r(n + l)r(-n~a) : 

o, 
1, 

0. 



if n e N, a £ Z_, 

if n G N, n + a G Z_ , 

if n G N, a G Z_, n + a G N , 
if n = 0, 
if n G Z_. 



(4.10) 



Note that for n G N, a G Z_ with n + a < -1, 



1 / d 
nl \ dx 



n + a 
n 

tl + OL 



-a- 1 



n + a 
n 



holds for n G No, Then Laguerre polynomials can be expressed as 

£=0 VJ 

for any a£l Remark that applying 14.11f> to the equation 

n-i 



By this definition, the equality 

(4.11) 

(4.12) 



1_ d_ 

nl \ dx 



1 



(n - 1)! 



1 / d 



! Vcfa; 



n+(a— 1) 
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with x = 1 and putting (3 = a + n, we have the identity 

fi ) = ( P - 1 ) + ( -\\ nEUEl (4.13) 
We introduce the polynomials 

Fjix) = -^Lf +1 {x), jeN , (4.14) 



Gj{x) = — \Lf +1 {x)- J -^—Lf_ x {x)}, j eN. (4.15) 



For k e N , j = 0,1,2, ...,k, let 



a k j = I ^ _ . I , it fc is even, 

fc + 2q /fc - 2 - j + b\ /fc-j + fA . 
ay = k _ 2 _. )-( k _. ), rffcuodd, (4.16) 



In Appendix IbI we will give the proof of the following lemmas. 
Lemma 4.1 For £ e No 



Lemma 4.2 For q,£ € 



with 



21 

F 2i (x)=Y / a njL v j (x), (4.17) 

2£+l 

G 2 £+i(.t) = ^ a 2t+hj L v j {x). (4.18) 



(F 29 ,G 2 , +1 )* = -<G 2m ,F 2? )» = r*<^, (4.19) 

(F2q,F 2 i)* = 0, (4.20) 

(Gaa+i.Ga/H-i), = 0, (4.21) 

_4r(2g + 2a + 2)_ 1) / 2g + 2a + lY ^ 



q (2« + l)l v 'V 2 9+! 

Then if we define the monic polynomials in a; of degree k for fc € No as 

Lemma l4~2l gives the following through the relation (|4.9|l and the orthogonality of the Laguerre polynomials 

JED. 



Lemma 4.3 For ?,<eN 

(-R 2g , R2£+l) = —(R2i+1, R 2 q) = r q S q £, 

(R 2q ,R2e) = 0, (R 2q +i, R21+1) = 0, 

where 

,-,ft\\ iq+1 (2g)!r(2 g + 2 + 2a) 



2 -2, T - K Ml ^ ^li^, (4.24) 

1 77 r(z/ + 1) 2 v ; 
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The choice of the polynomials Fj (x) and Gj (x) in (|4.14(1 and 14.15|l , and their explicit expansions in terms 
of the Laguerre polynomials fLemma |4.1|) are crucial, since they enable us to determine the appropriate skew- 
orthogonal polynomials ( Lemma 14. 3|) . As shown below, we are able to inverse the skew-symmetric matrix 
Aq given by (|4.1|l readily for arbitrary (even) TV, by using these skew-orthogonal polynomials. 

4.3 Matrix inversion 

Let &2fc = ^2fe+i = r k 1 ^ 2 , k £ No, and determine the polynomials {Mi(x)}o<i<N-i in (|3.12|) as 

M i (x) = b i R i (x), i = 0, 1, ■ • • ,N-1. 

Then by H4.1|) . I|4.4|l and Lemma f4. 31 we have the equality 

= (JN)ij, i,j = 1,2,- ■■ ,N, (4.25) 

where Jn — In/2 ® Ji- It is interesting to compare this result with l|3.7[) . Since J N = —In, we can 
immediately obtain the inversion matrix appearing in (|3.15|) as 

(i\(JN) A \j) = -(JN)ij, i,j = 1,2,- ■■ ,N. (4.26) 
If we consider a semi-infinite matrix 

J= lim J N =((i\J\j)) 

its inverse matrix may be given by 

j- 1 = (<*|i A ii>) f 

Using expansions lt5~6l and with Lemmas ETT1 IO and lUrl we can show 

(m,x\p\n,y) = (m,x\pJ\i)(i\J A \j)(j\n,y) 

and so 

(m,x\pJp\n,y) = (m,x\pJ\i)(i\J A \j)(j\Jp\n,y). 
Then the equations l|3.16[) are written as 

(m,x\pj\i)(i\(J N ) A \j)(j\n,y), if m > n, 

-(m,x\pJ\i)(i\(J A - (J N ) A )\j)(j\n,y) : if m < n, 
I m ' n (x,y) = (m,x\pJ\i)(i\(J A -(J N ) A )\j)(j\Jp\n,y). (4.27) 

Now we introduce the notations, just following the previous papers for multi-matrix models |33l 1151 I5T) . 



= -J. 



S m ' n (x,y) 



as 



i? 4 M (x) = L(m, x \i + 1) 

0% 



dy MyP^KtiMV^itrn-tuxly), (4-28) 



$| m) (x) = -~{m,x\pj\i + l) 

0* 



dyR < r\y)F^(y,x), (4.29) 
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for i = 0, l,--- , JV- 1,ot= 1,2,- •■ ,M + 1, where 
F {m \x lV ) = / diu / 



^«)(T-t m ,z|a;) ^ K )(T-i m ,u>|x) 
p(^ K )(T - t m) z|y) p(^ K )(T - t ro , to|yj 



(4.30) 



It should be noted that \x) = R i (x)p^ v ' K '>(t 1 ,x\0),0 < i < N - 1, and FW(x, y) = F(x,y), where 
Ri(x) and F(x,y) were defined by 14.23(1 and 14.3(1 . respectively. Then we arrive at the following explicit 
expressions for the elements of matrix kernel 13.18(1 of our Pfaffian processes, 



and 



(AT/2)-l 



£=0 

oc 



re 



1 



£=jV/2 1 



(JV/2)-l 



£=0 * 



(4.31) 



S™' n (z,y) = §%> n (x,y) = 5^™(x,2/) -^(i n -U,yk)l (m<n) . 



(4.32) 



5 Asymptotic Behavior of Correlation Functions 

In this section, we give the proof of our main theorem fTheorem l2.ll) . by estimating the N — > oo asymptotic 
of matrix kernel (|3.18|) of Theorem 13. II Elementary calculation needed for the estimation are summarized in 
Appendix IdI Here ajv ~ &at, N — > oo means ajv/bjv — * 1, N — > oo. We assume that T — N, t m = T + s m , 1 < 
m < M + 1 with si < s% < • • • < sm < sm+i = 0- We put 



Lj(x, -s m ) = L" (x) xin- and U<(x, s m ) = p^^^ j 0' 1 \ ,„' • ( : »- 1 I 

?( m ), 



5.1 Asymptotics of Rk(x) and Fc™'(x) 
Let 



Since c x ~ JV = T, 



i? 2 ,(x) - /(2€,b), 
JW(*) ~ _^_J(2^-l ) b)-7(2€ + l ) b-l)-J(2f,b-l) 

- |/(2£, b) - 2/(2£, b - 1), JV^oo, (5.2) 

where 

'(.■«)-g(*^r)^(5-)-ter)^(s-- 

for q G N and c S K. We set 

21 = NO, 
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and examine the asymptotic behavior of 1(21, c) as N — ► oo with some 9 £ (0, oo). When c e Z_, 



J 



can easily see 



-c- 1 

i 



Then from l|D.10jl in Lemma IDT2I with j = 2£ (i.e. r\ = 1 in l|D.5[) ). we 



I(2£,c) = J2 

j=0 



j + c 



(N0) 



c+u+l 



(e x ) 



—4- c - l \e,i,x,- Sl ), n- 



This result is generalized to the following lemma. 

Lemma 5.1 For any c E R, 8 € (0, oo), we have 

(N6) c+u+l 



1(21, c) 



Proof. 



21 

I(2£ >C ) = J2 

p=0 

2£ 

= E 

p=0 

21: 



p + c 
P 



^-p ( T7> — S l 



Jf-^^l.x.-si), N 



N 



p + c 
P 



p + c 
V 



{*Mj?.-0-*(£.-*)} + i: 

p=0 

t (T) (#.-) - 

p=0 \ * / k=0 p=Q \ V / 

p-1 1 /i\ 2f 

EEnr ^g,-, )+E 



2£ 

E 

p=0 



p + c 
P 



Repeating this procedure, we have 



p + c 
P 



^2£ T?) ~ S 1 



J(2*. c) = ^(-l) fe a fc (2^, c) Hl^, (|, -si) 



9 =0 



with 



a fe (2£,c) = ^ 

p=0 



P 



2/' 



2f / . x p-1 31-1 jfc-i-1 

^EE- E ' E 



ji=0j 2 =0 j fc =0 p=0 



p + c\ f p 



Using l|4.13|) . we can rewrite ak(2£, c) as 



It 



2£ 

E 

p=0 



P + c + 1\ {P + c\fp—l 

p - 1 



Of 



P-1 



2^ + c + 1\ /2£ 
2* 

Using this equation recursively, we obtain 

k 



o fc _i(2^-l,c + l). 



a A . 



(2£,c) = ^(-iy 



r=0 



21- 



1\ /2£ - 



2^- 



2£ + c + l\ (c + k 
21- k 



(5.3) 



(5.4) 



(5.5) 
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Thus JH3J with gi 



ives 



7(2^, C )=^(-1) 



'2€ + c + 1\ /c + // ' 



k=Q 



21 -k 



= Vy/ 



By simple calculation with the estimate (|D.2|) and (|D.10|) of Lemma ID. 21 we obtain (|5.3|l through the 
expression Ij2.16|l . | 

From above asymptotic of I(2£, c) with equations (|5.2() we have the following proposition. 
Proposition 5.2 (1) Suppose that £ e N and 2£ ~ iV#, AT — > oo for some 6 G (0, oo). T/ien, 



(N9) 2 a+1 



— 4-^(6, l,x,- Sl ), N 



oo. 



(2) Suppose that <eN and 2^ + 1 ~ N ^ oo for some 6 e (0, oo). Then 

aJi- b - l \e, 1, x, - Sl ) - J(-") (0, 1, - Sl ) 



2(iV6») [ '+ ,y 



We next examine asymptotic of R^ a \x). From the definition 1|4.28)1 and the expression l|4.7|) 

/* oo 

^i m) W=ci / driR k (c 1 T 1 p v ' K Ht 1 ,c 1 r,\Q)^ K \t m -t 1 ,x\c 1 r,) 



2"+ 1 i> + i) \xij \t 



ci\ ( 1 



a; exp 



' J J=0 



We put 



5(m) ^ _ 2"TT(i/ + l) f X i \ D M 



+ 1 + 2a) V ci 1 
If we set k ~ A^# as A" — > oo, (|D.1(I in Appendix [D] gives 



fcir 



i 



exp 



2T ~cZ, 



T(fc + l + 2a) V*r: 

then we obtain the following from Proposition 15. 21 

Proposition 5.3 (1) Suppose that i e N and 2£ — A^0, AT — > oo for some 6 G (0, oo). Then 

(2) Suppose that I G N and 21 + 1 ~ AT0, AT -> 00 for some 9 S (0, 00). T/ien 

aJ^-^ie, 1, x, -s m ) - J5- 6 ) (0, 1, x, -s m )l , Ar - 00. 



N 



(5.6) 
(5.7) 
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5.2 Asymptotics of 
Using q4.8p. (|4.30|) is rewritten as 

2(i/+l) 



oo oo / \ / 



Hence from (I4.29|) and l|5.6|l . we have 

/•oo 

^ m) (x) = c m / di/fljj™ 5 (^77)^(^77, a;) 



i) 



2"+ 1 r(i/ + l) Vxi/ \<W T« ^ 2Tc„ 

OO / \ 

where we have used the orthogonal relation IjB.lfl of Laguerre polynomials. Put 

* " ' = fc! I ~c~ x ) k " ' 

Then we have the following proposition. 

Proposition 5.4 (1) Suppose that £ G N and 2£ ~ NO, N — > oo /or some G (0, oo). Then 

/oo 
d?C°^ +1) («,^,Sm), N^OO. (5.10) 

(2) Suppose that I G N and 21 + 1 ~ NO, N —> oo for some 6 G (0, oo). Then 

. , 9/5—1+^^,^/2 

ft(x)- ^ ff+^fl.l.s.O, JV-oo. (5.11) 



Proof. 



00 



J<-i> x k/2 

3=0 p=0 

rp—U k/2 00 ^ 

3=0 ?=o 



Here we have used the notation IjC.ljl and introduced /3 = (J3j,k), which is the inverse of the matrix a = (ctkj) 
given by Lemma IC. II in Appendix IU1 By the skew orthogonality of {Qk} given by Lemma f4. 21 we have 

rp j, ft/2 * 00 

% } (*)~ \ r " E/w%(^*»)' * _f00 ' (5 - 12) 

*&(*) ~ - r "Y 2rl f>,^ (^,*») , iV->oo. (5.13) 

3=0 
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By jog and JCS, ^^Tl gives 

$S 1 (,)^-2r(2a + l)T-^( 2£ + 1 + + 1 2a ) 



oo 

X 



From ID. 81) we have 



A-2£-b-2W /a; N ~/ r -b-2W 



j=2£ v J 7 r=0 v 7 

r=0 ^ 7 p=0 



p J~2e+r+p \ N ^ m 



i — n \ / 



(5.15) 



where <|D.11|) of Lemma ID. 21 was applied. Setting j = 2lr\ — NO-q and using (|D.2J| in Appendix [D] we 
conclude that 



d? ? TThX^ ■ ^^oo- (5-16) 



-b-2+a 



NT(-b- 1 + a) ii (77-l) b + 2 -« 



Through the expression (|2.17(l . we obtain H5.11|l . 
By (jUl3)l of Lemma E3 with 



PjM+iLI (^,Sr. 



j=2f+l 

[tf+l)/2] 



j - 2r - b - 1 

6(1,2^+1) ^ ~'\N'" m J ^ ~'V j-2r + l 

; j=2£+l r=^+l v J 

1 sm, 



6(1,2*+ 1) 



where 



j — 2r — b — 1 
j - 2r + 1 



5(£)= £ 6(l,2r-l) ^(^' Sm ) 

r=t+l j=2r-l 

By this equation with the estimate l|D.3(l for 6(1, 2r — 1) and (|4.22() . I|5.12[l becomes 

~ -2T-(2^ + 2)-F(2a+l)^+ 1 + + 1 2a )W^W 

- -2T" l '(2* + 2)V J < /2 1 S'(£), TV -> oo. (5.17) 

From l(5~T5| with (|D~3ll 



r=C+l v 7 j=2r 



r^r ^fW) , 



2r(-b-l + a)y x A ' {r,-£,) b+2 - a 
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Thus 



fW 2 



i, ,a I j ^ (fit Vt x 1 s m) 



T{-b-l + a)J l ^ j t 1 



where we used the expression (I2.17[) . This completes the proof of Proposition [^] | 

5.3 Asymptotics of D m > n (x, y), I m ' n (x,y), S m ' n (x,y) and S m > n (x,y) 
From the expressions <|4.31ll with (|4.24() and the definitions (|5.7() and (|5.9() we have 

(N/2)~l ... 2o 



oo 
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-2a 



%' n {x, v )~- £ (i) *ST(*)«Sii(y)-*SSi(*)*S , (») 

e=(N/2) ^ ' 

(J\T/2)-l ^ 

, ' , ■■ X — x ~ 

=0 



1=0 

From Propositions 15.31 and 15. 41 we obtain the following asymptotics: 



D^' n (x,y) ~ V(s m ,x;s n ,y), 

fa ill) ~ 1( s mt x; s n , y), 
S^' n (x,y) ~ 5(s m ,x;s„,y), TV^oo, 

where T>,T,S are defined by l|2.19|l . 

Next we study the asymptotic behavior of p^ v ' K \t n — t m ,y\x). From 14.7J1 we have 



rt v > K Htn-t m ,y\x) 

v+l 



c rn a 1 I — j '/" >'X1> 



t/2 



t m X 

2Tc m 



exp 



-2 + ^)^ 
T I 2c„ 



Then by simple calculation with Lemma ID . 21 with a = 0, we have 

^ K) (t n -t m ,y\x) ~ ^) b/2 l^exp[2( Sm -s n )^] J v (2y/(fa)J v (2y/fajjj), 

X J=0 
/y\ b/2 

~ [-) 5(s ra ,i;S„,!/), TV ^ OO, 



where Q is defined by 12.21|l . and then S 7 ^' n (x,y) ~ S(s m ,x;s n ,y), N — > oo. Then, the proof of Theorem 
12. H is completed. 
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Appendices 



A Proof of (CD3D 



We assume that the number of particles N is even. Consider the multiple integral 

Af+l ,. M+l 



M 



n dx(m) det M J _ 1 (4 i) )^ K )(t 1 ,4 i) io)(i+xi(4 i) )) 

N(Af+i) i-J- Ki,7<JV 1 J ^ 

+ m=l 



X 



11 i j < w [^(Wi - i™,4 m+1) k (m) )(l + sgn (Mx (M+1) )) 



By the definition 1)3.4(1 with and jSH) , and by the equality ((3.11(1 . we have 



z n,tIx] 



[0] 



(A.l) 



where Zjy T [0] is obtained from Z^ T [x] by setting Xm( x ) = for all m = 1, 2, • • • , M + 1. 
By repeated applications of the Heine identity 



dx det 

l<i,j<N 



4>i{Xj) 



det 

l<i,j<N 



= det 

l<i,j<N 



dx 4>i(x)4>j (x) 



for square integrable continuous functions fa, fa, 1 < i < iV, we have 

M+l 



dy det 

l<i,j<N 
M ( 



/ TT d^" 1 ) J A/ l _ 1 ( 2 ;W)p(^)(i 1 ,a ; «|0)(l + xi^ (1) )) 



s + m= l 



Using the notations in Section it is expressed as 

ZnM = / rfy det (l + — |M+1,%) 

JR" 1<»J<W [ V 1 - XP+ J N 



dy det 

l<i,J<N 



(M + l, m| 1+px 



1 -p_x 



JV 



since (m, y) = (n, y\(p-) k \m, x) = for k > n — m > 0. Here we have used the Chapman- 

Kolmogorov equation, L dyp^ v ' K ^ (t — s, y\x)p^ u ' K ^ (u — t, z\y) = (u — s, z\x), < s < t < u < T, x,y & 
M+. Next we use the formula of de Bruijin [§] 



dy det 

l<i,j<N 



pf 



l<i,j<N 



dy / dysgn(y -y)fa(y)<f>j(y) 



for integrable continuous functions fa,l < i < N, in which the Pfafnan is defined by ((2.22(1 . Since (Pf {A)) 2 = 
det A for any even-dimensional skew-symmetric matrix A, we have 



z n,t[x] 



det 

l<i,j<N 



det 

l<i,j<N 



a i 



\ xp) J[l+PX-, ^—r) \j) 

1-XP+ J N V 1-P-X/jv 
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with 



(Ao)ij = (i|Jjv|i>, 

V 1-XP+ j n V 1-p+x 'N 



Since 



(A0« 

(As)ij 
= deti<j j< 



V-X> N 

(Ao)y , {PT> gives 



*Ar,r(f;0)^ = det 



% + (A) %)« + (A %)« + (A 



(A.2) 



By our notation l|3.13[) . (A )y = (i\(Jjsr) A \j), and it is easy to confirm that 1A.2|) is written in the form 



*£ r (f;0) 



det 

l<i,j<N 



(i\B\m,x)(m,x\C\j) 



(A.3) 



where we have introduced B and C as the following two-dimensional row and column vector-valued operators, 



B = 
C = 



pi 



The determinant (|A.3(1 is equivalent with the Fredholm determinant, 

Det(m,x|/ 2 +CoB|n, y). 

Introducing matrix- valued operators, 



we have 



1 - P+X 
1 

7 2 + CoB = K 1 
= K 1 



1 







1 -pJPx 



, K 



l-p-X 



A 



K_Kj 



pJJ N 



pj(l-j£j)p 



-J, 



N 



Jn Jp 



pJj£-p+ pJp-pJJmJp 



A 



-J, 



A 
JY 



Jm JP ~ P- 



K+ _1 K 



K+- 1 K, 



where J N = o( j A r) A o. From the orthogonality (|3.6J) and the definitions Ij3.9(l of the operators p+ and p-, 
we have the fact that 



Det (to, ar|K+|n, y) = Det (to, x|K_ |n, y) = Det (to, x[K\n,y) = 1. 



Then (|3.14|) is derived. 
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B Proofs of Lemmas 14.11 and 14.21 

We use the following orthogonal relations and formulae on Laguerre polynomials, which hold for a, f3 > — 1; 

T(a + 3 + 1) 



L?(x)LZ(x)x a e- x dx= >- 5 jk , j,keN , (B.l) 

x^L^(x) = jLJ(x) - (j + a)L?_ x {x), j e N, (B.2) 

^)=E( i_fc )^ a_1 )^ (x) ' jeNo - (R4) 

Remark 6. The identities ljB~2)> and ljB~3)l are given as Eqs. (6.2.6) and (6.2.7) in 0. The relation ijB~4|) 
is proved in [1] as (6.2.37) only when [3 > a > — 1. The identity (see (|^TU|l and g2|) 

A (I- a- \\ fk-t + a- 1\ . 

^{ i { k -e r m 

can be used to invert the relation l|B.4|) to the form 

3=0 V 

Therefore, the validity of (|B.4J| for /3 > a > — 1 implies that for a > (3 > — 1. 
B.l Proof of Lemma 14.11 

In this subsection we prove Lemma 14.11 which gives the expansion formulae of F k {x) and G k (x) in terms of 
{Lj(x)}. Taking the summation of the equalities (|B.3|I from to k £ No gives 

£ Ll\x) = -^Ltl^x) + ^Ll%x) = F k (x). (B.5) 

n=0 ' 

From ljB.4jl . (|B.5|) and l|4.13|l . we have 

F k {x) = L 2 k a {x)+F k _ 1 (x) 



k-j 



3=0 

fe /, . . „\ fe-1 



Since L%(x) = 1 and F (x) = Ll a {x) = 1, 

w - E ( fc fc i; b )w = w - E (* fc 1 7 ^ b )w 

= ^o(*)-(qW*) = 0. 
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Then we have (|4.17jl . From (|4.14jl and l|4.15|l we have 

G k (x) = -F k (x) + ^^F k ^(x), keN. (B.6) 
Then (|4.17(l gives (|4.18|l . This completes the proof. | 

B.2 Proof of Lemma [4.21 

We introduce a symmetric inner product 

/•OO 

(f,g) = / dxe- x x 2a f(x)g(x). 
Jo 

It is easy to see that it is related with the elementary skew-symmetric inner product (|4.5(l as 

/>oo POO 

(f,g)* = 2(If,g)- dxe- x ' 2 x a f{x) dy e-y/ 2 y a g(y), (B.7) 



where 

We consider the polynomials 



_ H dz e~ z / 2 z a f{z) 



? — x/2rg(l 



.?+2a r2a 



Wj[x) = Lf{x)- J —Lf_ 1 {x), jeN. (B.E 



Lemma B.l For j = N , k e N 



r( j + 2a + 2) 
(J + 1)! 



(W k ,F s ) = - ^ fe -i (B.9) 



Proof. For 2a > —1 the above equation can be derived immediately from the definitions (|B.8() . (|B.5ll . and 
the orthogonal relation (|B.1|I . We can extend it to the case a > —1. For k = 2,3, . . . it is easy to see that 



{Lf,x') = (-l) fe r(fc + 2o+l)^, j = l,2,...,fc. (B.10) 
However, if 2a < -1, (L|°, 1) = oo. Then we use the following equation: 

(Ll*,L?-L?(0)) = T{k + *« +1) 6 jk , keN , (B.ll) 



which is obtained from (|B.10(1 for j = 1, 2, . . . , A;, since L 2a (x) is a polynomial whose coefficient of the fe-th 
order is (— l) k /k\. By simple calculation we have 

(l,Lf -Lf(0))=0, ieN 0) 

and thus 

(L|<\ if - if (0)) = (L 2 k a if (0), Lf - if (0)) = 0, j > fc + 4. 
From the definitions of Wu and Fq = 1 it is easily to see that 

(W k ,F ) = -r(2a+2)4-i o. (B.12) 



2G 



When k,j E N, from ljR12)l and IRTTll . 



k\ 

This completes the proof. | 

Here we prove the following integral formula. 

Lemma B.2 For j e N, I e N 



(W k ,F 3 ) = Li* -t^Ll^L?) 

\ p=0 J 

p=l \ ' 
r(fc + 2a+l) 

- — <3fe-i j- (B.13) 



/ dxe- x/2 x a Gj(x) = 2e _1!/ VW , J -(a:), 
Jo 

/ dxe- x/2 x a F 2e (x) = 2 a 
Jo 



7 (a+l,z/2) 



r 

r=0 



where 7(0, y), c > is the incomplete gamma function 7(0, y) = dxe x x c 1 . 
Remark 7. If we set a = in (|B.14J| . we will have the simpler equation 

[ dx e ~ x/2 Tx{ L ^ {x) L °^ {x) ) = 2 ^ Z/2 [^M) ~ Q-ii*) 



J'GN. 



Proof of Lemma iB.Sl We first introduce the functions defined by 

1#°(*0 = £ dxe-^ 2 x a ^-Lf(x), 36% 
</> 2o (z) = f dxe~ x ' 2 x a Lf (x), j G No. 



(B.14) 



Then (jB73|) gives 

^00 = "^*) +$"(*)■ (B.16) 

On the other hand, by l|B.2|l . 



= j / dxe-z^x^Wjix). 
Jo 

Noting the assumption a > — 1 and the fact that Wj(x) = 0(x), in x — > 0, 

2a^f (z) - 2j je-^VW^z) + jf da;a; a Qe"^ 2 ^^) - e ~ X/2 ^ W i( 3 

= 2 3 e' z ' 2 z a W 3 (z) + jyfiz) - (j + 2a)yf_ x {z) - 2^f (z) + 2(j + 2a)$'iC*). 
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Then we use (JB.16I) to eliminate <p 2a (z) and (p^_ x (z), and the equality 1B.14|) is obtained from the relation 



dx e- x ' 2 x a G 3 {x) = 4>f +1 {z) - 3 -^^f_ 1 {z). 
o J 



Note that (|B.14J| gives a recurrence relation for ip 2a (z), 



tfli(z) = ^T^i-i(*) + 2e-^z"W J (z). 



It is solved as 



t + a\ iA ft - r + a\ 1 



W 2 i-2r(z). 



^+i(*) = (^ a )v^W + 2e-/v(^ a )x: 
Since L 2a (x) = (1 + 2o) - x, 

iPl a {z) = [ dxe- x/2 x a ^L\ a {z) = - ( dxe~ x/2 x a = -2 0+1 7 (a + 1, z/2). 
Jo dx Jo 



Then we have (|B.15I) . | 

Then we can prove Lemma l4~2l From Lemma fB . 21 and two relations (|4.15(l and (|B.7f) we have 

(F2 q ,G2i+l)* = —{G2l+l,F2q)*——^{W2i+l 1 F2 q ) 1 

(G 2q +i,G2e+i)* = 4(T / t^2<j+i , G21+1) 



Then (|4.2U|) and l|4.21|l are derived from Lemma fB. II From Lemma IB . 21 and l|B.9 

X) 

dw e~ w/2 w a j{a + 1, w/2)F 2i {w) 

7 (a + l,w/2) / dze~ z/2 z a F 2e (z) 
Jo 

dw 2- ( - 1+ ^e- w / 2 w a / dze~ z l 2 z a F 2 t{z) 
Jo Jo 

r+1 C + i °) { r(a + ~ So dW 2 ~ {1+a)e ~ W/2wa ^ a + w /2) 



2 a+l 



£ + a 



( £ + a ) {r(a + I) 2 - \ [ 7 (a + 1, w/2)X } =^Y) n « i 1 r 



Then by using Lemma IB. 21 and l|B.9(l we have H4.20|) , since 



(F 2q ,F,. 



2< I 



dwe' w/2 w a F 2 t{w) / dze- z/2 z a F2q{z) 



dwe- w/2 w a F 2e (w) / dze- z l 2 z a F 2 q{z) 



1 J Jo 
,2a+2 /^7 + a 
9 



e- lo/2 w a 7(a + 1, w/2)F 2e (w) 



r(o+ l) 2 = 0. 



This completes the proof. 
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C Inverse of {ctk,j} 

We put 

Q2t(x) = F 2 i{x) and Qu+i(x) = G 2 e+i(x), (C.l) 

for I 6 No, and Qk = for fe 6 Z_. Lemma f4. II gives the expansion formula of Qk(x) in terms of {Lj(x)}. 
Here we give the formula to expand Lj(x) in terms of {Qk{x)}- In other words, we calculate the inverse of 
the matrix a = (ctfc.i) given by i|4.16[l . which is denoted by (3 = (flj,k) and used in Section 5.2. 
Let 6(n) = (n + 2a) /n, neN, and 

!b(m)b(m + 2) • • ■ 6(n), if to, n are odd and m < n, 

1, if TO,rt are odd and to > n, (C.2) 
0, otherwise. 

Then the following lemma holds. 
Lemma C.l For j 6 No 

L»{x) = Y,p3,kQk{x). (C.3) 

fe=0 

where (3j,k, k,l,. . . ,j, j € No are defined by the following: 
When k is even 

( 0, ifj < k, 



and, when k is odd 



0, ifj < k, 

- 2r - b - 
j - 2r + 1 



I [0'+l)/2] / ' •> r I X 

^ = < - E 6(* + 2,ar-i)f J T^7 1 ), ^i>*. (c ' 5) 



r=[(fc+l)/2] 
Proof. From l|B.5(l and ()B.6|) we have 

11 2£+l 21-1 

Q 2l (x) =Y,Lf{x), Q 2 i+i(x) = -J2 Lf{x)+b(2l+l) E if(aj). 
3=0 i=o j=0 

By simple calculations we see that 

fc=0 

where /3 - fc , k = 0, 1, . . . ,j, j € No, are defined by the following: 
When fc is even 

1, if j = k, 
P j ,k={ -1 ifj = * + l, 



and, when fc is odd 
Using the formula l|B.4|l . 1C.6(I gives 



0, otherwise, 

b(k + 2,j — 1), if j is even, 
^ fe = \ -b(k + 2,j), ifj is odd. 
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If we define 

3 

E' 

p—k 

for j > k and /Sj^ = for j < k, l|C.3j) is satisfied. The expressions (IC.4|l and l)C.5(l are derived from (|C.7(I 
by simple calculation. | 



D Elementary Calculation for Asymptotics Estimation 

By Stirling's formula T(x) ~ v2nx x ~ 1 ' 2 e~ x , a; — > oo, we have 

F(n + a + l) 



r(n + 1) 

/ n + a\ (n + 1) 



n / r(a+l) 



+ ji^oo, (D.I) 

-> oo, (D.2) 



for any a £ R \ Z_ , and 



2 a r(a + i)V 2i 

^ + l,2p-l) = |£^~(!V, /^oo (D.4) 



6(1,2^-1) \l 

for £,p with £ < p, where &(m, n) is defined by (|C.2J) . 

From now on, we assume that T = N, t m =T + s m with s m < 0. We set 

2£ = NO and j = 2£r), (D.5) 

and consider the limit N — > oo with some r],0 G (0, oo). Then we have 

9T-f V / 9« \ ^ 



Xm = ( — 7 — — J = ( 1 _ J ~exp(-2fl m ^), AT^oo, 
and 

^ ( " 1)P (p) x ™ P ^ (2£ra (i) Qexp( " 2Sm ^ ) ' N ^°°' (D ' 6) 

We use the following identities (see (|4.1U|) and pages 8, 201 and 202 in |32*]1. 

(i) Let a G N and c e R. Then 

±(-ir( a )( n - p+c ) = ( n - a+c ). (d.7) 

^ \Pj\n-p-jJ \ n-j J 

(ii) Let a e No, c £ K and a^, fc = 1, 2, . . . , be a sequence in R. Then 

£(T)*-i:(' + r +a )£<-iK:W (D ' 8) 



(iii) Let a 6 No, and a/., fc = 1, 2, . . . , be sequences in R. Then 

, a \ a-fi 



fc=0 v 7 /3=0 ' p=0 ^ / q=0 
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Lemma D.l For any a S No and w > we have 

Proof. From the definition of the Laguerre polynomials (|4.12|l and l|D.7jl . we have 

n — p + v \ y J 
. n — v — j ' ) j ! 



x><K» - g(- 1 )<;)g ( - 1 y(:::!^ 

(— 1) J fn — a + iA 



E 

J=0 



j! \ n-j 



y 3 ■ 



Hence, by CE2 



oo 

£(-*) 



«A .. \ ^ (— 1) J /n — a + zA fw\i 



n— >oo \ n / 1 — ' ' W J " '' \ n J n^oo * — ' 7! \ n — 7 / \n 

— a. 



i+V — OL 



1. w 



T{j-a + v+l)j\ 



dw ' \fc5 r(i/+i + i)j! 

Then we obtain the lemma | 

Applying the above lemma, we obtain the following asymptotics, where and are defined by (|5.1(l . 
Lemma D.2 For any a € No, n € (0, oo), and a; G K+ we have 

E(-l) P ( a )£- +P (^,s m ) ~(2£)- a jV>(9,ri, x,s m ), iV^oo. (D.ll) 



p=0 



Proof. Since 



A.m 



±G)i><%(*)lH a .O 



p=0 

a-/3 



Y 3-q 

Am 7 



by (|D.9|I . the asymptotic (|D.10|) is derived from (|D.6|I and Lemma fD. II with n = j = N6r),w = dr/x. From 
HD.9fl . we have 



fc=0 v 7 



(D.12) 
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By (|D.1|) . we see 

a-0 / o\ w • ,< a— 0—1 



^ Us ^(i + g + l + z,) " ^ 1 U 9 )r(j+ q + 2 + v) 

v{v + \)---{v + a- p-l)T{] + \) 
T(j + 1 + a - + v) 
/ A \ a ~ (i 

~ ( 2 ^)-(«-/3+-) _A jV->oo. 



f/// 



On the other hand, (|D . 1 0|) gives 



^ ( " 1)P (p)^ +P (^) Xm ° +P) ~ (2 ' r_/3 ("^) K^^} 



Hence, the asymptotic IjD.llll is derived from ID. 12(1 . | 

E On Temporally Homogeneous Limit 

Lemma E.l For any c G R and tj,9,x > 0, we have that as t — > oo 

J^{6,r],x,t) ~ {2t9)%6r]x) v ' 2 J v {2^/6rix)e' lt07 ', (E.l) 
J^ c) (6»,?7,a;,-t) - (2t9) c (6r]x)-^ 2 J v {2 V^K 2 * ", (E.2) 

d£ C a ^ c+1) (^,^^ ; -i) ~ (2^) c (^)-" /2 J„(2V^)e- 2te . (E.3) 



l 



Proof. From the expression (|2.16|) with the definition H2.12fl . we have 



4 C) (M,M) = ^Eli^-E (J)^(M,*,0)(2*)*-J 



{2teyj v {e, v , x , o) E ( ffi^ ■ 



r(-c)^ /-v— - ^ fc!(fc _ c) 



From the relation 



and the equation 



(see (1.1.19) in 0), we have 



E 



d ^ (-l) fc z fe - c 



e?z ^— ' k\(k — c) 



00 / i\fe roo 



-<->e- 



r( _ c)=lim f;(=Jffl£: 

t^oo^ fc!(fc-c) 



/c=0 

Then we conclude 

Ji c \e^x,t) ~ e 2te "(2^) c X(^,r/,a;,0) = (2te) c J u {6,r],x,t), t 
Hence (fE.lfl is derived from (|2.12(1 . 
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Let n = [c + and (3 > with c = n — [3. Since 

X(e,T],x,-t) ~ (2t6) n J„(6,r],x,-t), t^oo, 



d 
dr) 



fTTTjl gives 



V2m "~' r dCC f3 - 1 l(0,v + ^^-t) 



Then (|K2| is derived from lj2~T3|) . 
From (IE.2II we have 



r(/3) 7 * * v ' ' 2*0 
(2t6) c J v (6,r),x,-t), t^oo. 



/OO /-GO 

- (2t0) c ^,(0,l,x,-t), t->oo. 

This completes the proof. | 

Applying the above lemma, we have as s m , s n — > — oo with the difference s n — s m fixed 

(zy) b / 2 (s„-s m ) /- 1 



(sm + SnX^Z/) 6 / 2 
and 

6/2 /- 1 



S(s m ,x;s n ,y)~(T) b/2 d0 J„{2^)J u {2^Wy)e 



2(s m -s„)6 



It is then clear that 



lim P(a m , a;; s„, y)l(s m , x; s n , y) = 0. 
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